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Abstract 

In this paper, we consider generic coranlc 2 sub-Riemannian structures, and we show that 
the Spherical Hausdorf measure is always a C^-smooth volume, which is in fact generically C^- 
smooth out of a stratified subset of codimension 7. In particular, for rank 4, it is generically 
C^. This is the continuation of a previous work by the auhors. 

subjclass: 53C17, 49J15, 58C35 

Keywords: optimal control, sub-Riemannian geometry 

1 Introduction 

In this paper we consider sub-Riemannian structures s = (A,g) over an oriented n-dimcnsional 
manifold M. The distribution A has rank p and corank k = n — p, and g is a Riemannian metric 
over A. In most of the paper, k = 2. Moreover, the distribution is assumed to be 2-step bracket 
generating. 

The set S of such (corank 2, 2-step bracket generating) sub-Riemannian structures over M is 
endowed with the C°° Whitney topology. 

As it will be recalled in the next section, there is a natural smooth measure associated with 
the structure s, called the Popp measure (see [H])- It has been shown in [T] that the Radon- 
Nykodim derivative fsviO of the spherical Haussdorf measure with respect to the Popp measure 
at a point ^ is just (universally) proportional to the inverse of the Popp-volume of the unit ball 
of the nilpotent approximation of s at ^. Moreover, in the same paper, when fc = 1, it is shown 
that fsv{£,) is a function {C^ along curves), which is not in general. 

The nonsmoothness appears only in sub-Riemannian structures for which the nilpotent ap- 
proximation depends on the point and can occur at points (called resonance points) where certain 
invariants of the structures coincide. The (high) degree 3 of differentiability is due to the fact 
that, in the corank 1 case, the conjugate locus of the nilpotent approximation coincides with the 
cut locus. This coincidence is no more true for higher corank. In particular, this is shown in [4], 
in the corank 2 case, and an explicit characterization of the cut-locus is given. In the same paper, 
as a simple byproduct, it has been shown that fsv is generically for p = 4, fc = 2. 
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Starting from the explicit characterization of the cut-locus obtained in [J , in the current paper 
we go one step further and obtain the following result. 

Theorem 1 (step 2, corank 2 ) We have the following: 

1. the Radon- Nykodim derivative fgv is always ; 

2. The Radon-Nykodim derivative fg-p is genericall^ , out of a stratified set of codimension 
7. In the particular case p — 4, there is an open-dense subset of S for which fsv is -smooth. 

Remark 1 In the case of a non-orientable manifold the Popp measure cannot be defined as a 
volume form, but just as a density. However, Theorem [I] still holds true since it is essentially 
local. 

Roughly speaking, fs-p depends on the maximum eigenvalue of a certain skew symmetric ma- 
trix (depending on the point) defining the nilpotent approximation of the structure at the point. 
(This eigenvalue is an invariant of the structure.) Hence, the study of the differentiability proper- 
ties of fsv requires a fine analysis of the regularity of the maximum eigenvalue of a family of skew 
symmetric matrices smoothly depending on parameters. When the maximum eigenvalue is simple 
at a point, then in a neighborhood of that point it is C°°. A drop of regularity appears at points 
where the maximum eigenvalue is multiple. The regularity can be obtained as a consequence of 
the fact that when the maximum eigenvalue is multiple, the cut locus coincides with the conjugate 
locus. H This fact does not permit to get the result, which requires a deeper analysis. To treat 
double eigenvalues we need an adaptation of a deep result of Arnold [3] , to the case of versal defor- 
mations of real skew-symmetric matrices. The case of triple eigenvalues is apparently extremely 
difficult and we do not treat it in this paper. However, the set of skew symmetric matrices with a 
triple eigenvalue is an algebraic subset of codimension 8 in skew symmetric matrices (we provide 
a proof of this technical fact in appendix) . In the particular case of rank 4 this set is generically 
empty. 

The paper is organized as follows: in Section [2l we recall the definition of the Popp measure 
and that of the nilpotent approximation of s = (A, g). We avoid to recall all standard definitions of 
sub-Riemannian geometry since these are already given in [4]. Then, we recall the main result of 
[3] which gives the cut time for geodesies issued from the origin. This is our key point. In Section 
|3l we give the proof of Theorem [TJ This proof uses a certain number of technical tools that are 
collected in appendix. In Appendix I A. 1 1 we recall certain basic facts about quaternions, which here 
represents a very convenient tool. In IA.2I we study versal deformations of real skew-symmetric 
matrices. In lA. 31 we discuss the codimension of the set of skew symmetric matrices having either 
a double or a triple eigenvalue. In IA.4I we prove a result that (generically) allow us to make a 
crucial change of coordinates. In lA.51 we recall how to get a useful formula for the volume of the 
nilpotent ball in the corank 2 case. 

2 Prerequisites 

2.1 Nilpotent approximation 

-'In this theorem, genericity means that the property is satisfied for a subset of sub-Riemannian metrics that 
CONTAINS an open-dense set. Indeed, in the transversality arguments, we can always avoid the closure of certain 
Whitney-stratified bad-sets, in place of avoiding the bad-sets themselves. 

■^Notice that in the corank 1 case, the cut locus coincides with the conjugate locus at every point, see [l]. 
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We define the nilpotent approximation in the two-step baracket generating case only. The tensor 
mapping: 

[.,.]: A^xA^^T^M/A^, (1) 
is skew symmetric. Then, for any Z* e {T^M / A^)* we have: 

Z*{[X, Y] + A^) =< Az^ (X), Y >g 

for some (^-skew-symmetric endomorphism Az* of A^. The mapping Z* — > Az* is hnear, and its 
image is denoted by C^. 

The space = A^ T^M/A^ is endowed with the structure of a 2-step nilpotent Lie algebra 
with the bracket: 

m,Wi), iV2,W2)] = (0, [Vi,V2] + A,). 

The associated simply connected nilpotent Lie group is denoted by Gj, and the exponential map- 
ping Exp : — >■ is one-to-one and onto. By translation, the metric over A^ allows to define 
a left-invariant sub-Riemannian structure over G^, called the nilpotent approximation of {A,g) at 

Any fc-dimensional vector subspace of T^M , transversal to A^, allows to identify and G^ 
to T^M ~ Aj © TjM/Aj. If we fix G M,we can chose linear coordinates x in A^j, such that the 
metric g^^ is the standard Euclidean metric, and for any linear coordinate system y in V^q, there 
are skew-symmetric matrices Li, G so(p, R) such that te mapping [T] writes: 



[X,Y]+A^ 



( X'LiY \ 



V X'LkY J 



where X' denotes the transpose of the vector X. 

This construction works for any A, but A is 2-step bracket generating iff the endomorphisms 
of A^, Li, i = 1, fc (respectively the matices Li if coordinates y in are chosen) are linearly 
independant. 



2.2 Popp Measure 

In the 2-step bracket generating case, the linear coordinates y in T^M/ A^ can be chosen in such a 
way that the endomorphisms L^, i = 1, fc are orthonormal with respect to the Hilbert-Schmidt 
scalar product < Li,Lj >= -Tra,ceg{L'^Lj). This choice defines a canonical euclidean structure 
over T^M/A^ and a corresponding volume in this space. Then using the Euclidean structure over 
Aj, we get a canonical eucildean structure over A^ © T^M/A^. The choice of the subspace 
induces an euclidean structure on T^M that depends on the choice of V^, but the associated 
volume over T^M is independant of this choice. 

Definition 1 This volume form on M is called the Popp measure. 

By construction, the Popp measure is a smooth volume form. 

Let us recall a main result from [T]. 

Theorem 2 (equiregular, any step, any corank) The value fsviO o,t(^ Cz M of the Radon-Nykodim 
derivative of the spherical Hausdorf measure with respect to the Popp measure is equal to 2^/fl{B^), 
where Q is the Hausdorff dimension of the sub-Riemannian structure as metric space and fJ-iB^) 
is the Popp volume of the unit ball of the nilpotent approximation at ^. 



3 



2.3 Geodesies and Cut-locus 



We restrict to the corank 2 case. Here, we consider geodesies of the nilpotent approximation of 
s = (A,g) in T^gM ~ R", issued from the origin. A transversal subspace V^^ is chosen, together 
with the linear Hilbert-Schmidt-orthonormal coordinates y in Vj^, and euclidean coordinates x in 
A^,,. The geodesies are projections on R" of trajectories of the smooth Hamiltonian H on r*R" : 

p 

H{p'',py,x,y) = sup + Vpfwj +pfa;'LiM + pfa;'i2u)- (2) 

where are the coordinates dual to x, y. Geodesies are arclength-parametrized as soon as the 

initial covector {p^ (0) , (0)) verifies 7J(p^(0),p^(0), a;(0), j/(0)) = 1/2. For geodesies issued from 
the origin, this condition reads ||u(0)|i = |b'^(0)|| = 1, where the norm || • || is the one induced by 
duality on A|^^ . 

Note that p\,p^ are constant along geodesies, since the Hamiltonian ^ does not depend on 
the ^-coordinates. 

The following result is shown in [4], and is crucial for the proof of our result. 

Theorem 3 The cut time tcut of the arclength-parametrized geodesic issued from the origin and 
corresponding to the initial covectors (p^(0),p^) is given by: 

tc'i 



max(cr(p^ii +^^^2)' 

where max((T(A)) denotes the maximum modulus of the eigenvalues of the skew symmetric matrix 
A. In general, the conjugate time is not equal to the cut time. 

Remark 2 In fact the cut time is also conjugate if and only if the matrix p\Li + pf L2 has a 
double maximum eigenvalue or [Li,L2\ ~ 0. 

It turns out that the singularities of the Hausdorf measure appear due to collision between the 
two largest moduli of eigenvalues of the matrix p\Li +^2-^2- The set of skew-symmetric matrices 
that have a double eigenvalue is a codimension 3 algebraic subset of so(p, M) (see Appendix IA.3|) . 
Then, from the tranversality theorems ([2]), for generic (open, dense) sub-Riemannian structures, 
the set S2 of points of M such that PiLi +P2L2 has a double (at least) eigenvalue for some Pi,P2 
has codimension 2 in M. The problems of smoothness of the Hausdorf measure will occur on E2 
only. 

Along the paper we set, for the geodesic under consideration: 

2tt 

p\ = rcos{0), pI = rsin(6'), ^5(6*, r) = rm — \ — yT~v^^ wheie £, = {x,y) G M. 

It is known ([51[71ll]) that A^{9,r) is a Lipschitz function of all parameters £,,0,r. We write 
also A^{9) = A^{0,1). 

3 Proof of Theorem [1] 

For a fixed point ~ {xo,yo) G M, let us consider the exponential mapping £ associated with 
the nilpotent approximation at where x,y are coordinates as in Section \2l2\ 

f,(pg,Po^)=^(e*^(pg,p^,eo)), 
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where tt : T*M — > M is the canonical projection, and H is the Hamiltonian vector field associated 
with the Hamiltonian ([2]). Here {pq,Pq) arc initial covcctors satisfying H{pq,Pq,^o) = 1/2. 

As above we have (i) = Po = (PijPI) = {r cos{6) , r sm{6)) . Also, by homogeneity, £t(pg,pQ) = 
£i{t pl.tpl). 

In our paper [4] , the following formula is given for the volume at a point ^ e M of the unit 
ball of the nilpotent approximation. For the benefit of the reader, this formula is established here 
in Appendix lA. 51 

n - / / / drde (3) 



^0 



where B is the unit ball in the euclidean p-dimcnsional pg-space, and J{pq, 0, r, ^) is the jacobian 
determinant of £i(pg, r cos(6'), r sin(6')). 

We set /^(6i, r) = J(p§, 61, r, O^pg, and W^i0) = /^{O, r)dr. If wc show that 1^5(6') is 

or w.r.t (e,^), it will imply that is or w.r.t ^. 
In a neighborhood of a fixed (0o,Co) & y. M we have, 

W^{e)^ / U{e,r)dr (4) 
Jo 

^50 (So) Ms (9) 

/^(0,r)dr + / hiO,r)dr 

(/) + {II)- 

The term (I) is smooth. We are then left to study the smoothness of II{^,9). 
3.1 Proof of the fact that W^{e) is always 

Setting z ~ zq = {Oo,(,o), and f{z,r) ~ f^{9,r), A{z) = the tangent mapping to 

1 1(^,0), at {00, ^0) is 

D II{zo){h) = Y^f{z,,A{zo))^{z,)h.,. (5) 



i=l 



This last expression makes sense, and is continuous w.r.t zq for the following reasons: first as 
we said, A{z) is Lipschitz-continuous, then the derivatives are bounded. Moreover at points 
such that A is not differentiable, /(zq, A(zo)) vanishes. This last point follows from the fact that 
when the eigenvalue of A{zq) having maximum modulus is multiple then the conjugate time is 
equal to the cut time, which makes the jacobian determinant J{pq,0q, A{0o, Co), Co) vanish for all 
Pq. This comes from the section II. 3 1 in the paper [T]. 



Remark 3 In fact, it follows from the same paper that, if A{zq) corresponds to a multiple eigen- 
value, then the rank of J^^{pq,0(),A{0q,^q)) = J{pq,0q,A{0o,£,o)-,Co) drops by 2 at least, indepen- 
dently of Pq. This point will be very important in the next section. 

This ends the proof. 
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3.2 Proof of the result 



It follows from the transvcrsality theorems ([21 [5]) and from Lemma [5] and Lemma [3] in the 
Appendix, that there exists an open dense subset of sub-Ricmannian metrics, still denoted by S, 
such that all elements s of 5 meet: the set Us C x M of (6*,^) such that A^{d) corresponds 
to a triple (at least) eigenvalue is a locally finite union of manifolds, regularly embedded, of 
codimension 8 in S"^ x M, and the set Us C x M of {0,£,) such that A{6,^) corresponds to a 
double (and not triple) eigenvalue is a locally finite union of manifolds, of codimension 3. 

We want to show the following property (P), for a (smaller) generic (residual in the Whitney 
topology) set Sq of sub-Riemannian metrics over M: 

(P) the partial derivatives Di{z) = f{z, A{z))^{z) from ^ are in a neighborhood of all 
points zq such that A{zq) corresponds to a double (and not triple) eigenvalue. 

To do this, we fix sq and zq € Usg and we consider a (mini)versal deformation of L(^Oj^o) = 
ii(Co) cos(^^o) + -^2(^0) sin(0o) = L{zo), as introduced in Appendix lA. 21 It follows that: 

Li(,e)^Liz)^giz)-'Tiz)giz) 

where g{z) belongs to the orthogonal group an^l the functions g(-), A(-), g(-), A(-) are smooth with 
respect to z. 

The following crucial Lemma is proved in Appendix IA.4I 
Lemma 1 The property 

(R): the map x M 3 z t-^ q{z) G M.^ , has rank 3 at every z £ Us, 
is residual in S. 

Let us call Sq the subset of S for which (R) holds. If sq is fixed in Sq and Zq G Usq then, 
locally around zq, we can find a system of coordinates in 5^ x M in such a way that the three first 
coordinates, 2:1,^2, 23 become the three components of q{z). Note that these 3 coordinates vanish 

at Zq. 

Locally, the codimension 3 manifold Us„ is determined by the equations zi = Z2 = ^3 = 0. 

As we said in Remark|3l the rank of J(p§, z, A{z)) drops by 2 at least, indepcndantly of pg, at 
each point z € Usg- Formula ([7]) in the appendix tell us that A{z) = -S= — | — = where A(z) 

is smooth and nonzero. We set Z4 = (24, Zn+i) and zi — (zi, Z2, z^). 
The Jacobian determinant J(pg,z,r) can be written as 

Vi{po,z,r)A...AVn+i{po,z,r), 

for certain smooth n + 1-dimensional vectors Vi(pq, z, r) . . . V^„+i(pg, z, r). 

For all Pq, at points (z, r) such that zi — 0, r = A{z) = the vectors Vi . . . Vn+i have rank 
?i — 1 at most. Then 

'^T(z) plays the role of T{fi(z)) in Appendix IA.2I Tiz) is the block-diagonal matrix T{z) = Bd(A(2:)g + 
q{z), A{z)). Here, following the notation introduced in the appendix, g is a pure quaternion, 5 is a pure skew- 
quaternion, A{^) is a 2 X 2 block diagonal skew-symmetric matrix and A{^) is a nonzero real number. 
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dJ dV dV 
Q^T^ = ^hV2h...h Vn+i + l/i A — ^ A . . . A 14+1 + . . . 
ozi ozi ozi 

and 

dJ dVi dV2 

0= — = —AV2/\.../\ Vn+l + l/l A — A . . . A Vn+l +■■■ 

or or Or 
It follows that J, 1^, 1^ vanish at all (z, r) with li = 0, r = A{z) = 

Therefore f{z,r) = Jgi J^{pQ,9,r) dpg is a quadratic expression in the variable ii, r — -^|y 
depending smoothly on z, r: 

f{z,r) = Q,.r{zi,r - j^). (6) 

Now we study the continuity of the second partial derivatives of W^{0) — f^(^0,r)dr 

+ ZaSSo) '^)^'^' ^it^ notations, W^(z) = f{z, r)dr + /^/^^^ /(z, r)dr. 

The first partial derivatives, at any point zq were: 

dW ^ M(-o) 3 «^ 

^'^0)" / -^l[zo,r)dr + i{zQ,A{zQ)) — {zQ), 
OZi Jq OZi OZi 

=:III{zo) + IV{zo) 

To show that ^^glj^'^ exists and is continuous, we proceed exactly as in Section |3.1[ using the 
fact thatT^/(z,r) also vanishes at (zi = 0, r = -rri)- 



The more difficult point is to show that ^"^^^.^^ exists and is continuous. 



(£ 

dzj 



dzj dzj \ ' dzi 



We get: 



^^^-(z).|^(z,A(z))^)4-?^(z,A(z))^^4-/(z,A(z))-^^^(^) 



dzj dzj dzi dr dzi dzj dzidzj 

=: V{z) + VI{z) + VII{z). 



The casesof F(z), F/(z) are obvious, since again is bounded, and the functions §j-{z, A{z)), §^(z, A{z)) 
are continuous and go to zero when zi tends to zero. The only difhculty is the case of VII{z). 

Remind that A[z) = ^^^^^Jj -^|| where A(z) is nonzero, smooth. Then the only problem may 
occur for z = 1, 2, 3. 

Let us consider only the 2 cases: (1) i = 1, j = A, (2) i = l,j = 2, the other being similar. 
Case (1): ^ = (,(,)+f,,||). (^+p^), and £fg is bounded. It is multiplied by /(z, A(z)), 
which tends to zero when zi tends to zero. Then it is zero at points zi = 0, and it is continuous. 
Casc(2):^^ = _z2vr / sa ^ ^ _^ £1(^)^122 ^herc C(z) is 

^ ' ozi (A(z) + ||zi ^ 021 Ipill'' dzidz2 \ J ' \ v / 

bounded, D{z) is continuous. Then, the question is the continuity to zero of (p{z) = ^^^j^^jp^, in a 
neighborhood of the set E ~ {zi = 0}. Let us use Formula It gives /(z, A{z)) = Qz^r{zi, A{z) — 
-^^)- But A{z) = x{z)+\\zi\\ ' then, A{z) — = ^{z)\\zi\\, where ip{z) is continuous. It follows 
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that ip{z) tends to zero when zi tends to zero. The sub-Riemannian volume is in a neighborhood 

It follows that fsviO is generically except on a bad set of codimension 8 in x M, and 
the theorem is proved. In the case n = 6, the bad set is generically empty in x M and property 
(R) is open dense in S. 



A Appendix 

A.l Pure Quaternions in so(4) 

In so(4), it is natural and useful for computations to use quaternionic notations. Set: 
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The matrices (resp. i,j, k) generate the so-called pure quaternions (resp. pure skew- 

quaternions), the space of which is denoted by Q (resp. Q). The Lie algebra so(4) = Q ® Q, and 
quaternions commute with skew-quaternions: [Q, Q] — 0. 

We endow so(4) with the Hilbcrt-Schmidt scalar product: < Li,L2 >= trace{L[L2)- 
Then, k, k form an orthonormal basis. The eigenvalues wi, ^2 of A = q + q meet: 

-M' = i\\q\\±\\q\\?- (7) 
As a consequence, an element A E so(4) has a double eigenvalue iS A E QU Q. 

A. 2 Versal deformation of skew-symmetric matrices 

The results of Arnold in [3] can be easily extended to the real smooth case {C°°), for skew- 
symmetric matrices, under the action of the orthogonal group: 

Theorem 4 Let N{p) be a family ofnxn matrices smoothly depending onp at (R',0). Let On 
be the orbit of N = N{0) under the action of Gl{n,W) by conjugation. LetT{pL) be a smooth family 
of matrices, depending on the parameter jj. £ M.'', such that the mapping ji — > T(/i) transversally 
intersects On at some N = g~^Ng. Then, there is a family of (smoothly depending onp) matrices 
g(jp) and a smooth mapping p — > such that N{p) = g{p)^^T{^(ji))g(p). Moreover, for the 

transversal T{p), one can choose the centralizer of N in gl{n,M.). 

We rephrase the result in the case of a skew-symmetric matrix N that has a double (but not 
triple) eigenvalue. Then, by section ETTl we can assume that N is (conjugate to) a block-diagonal 
Bd{aq,S), where q is a unit skew-quaternion and 5 is a block-diagonal skew symmetric matrix 
with 2x2 blocks and non multiple eigenvalues. The centralizer of q in so(4, M) is the vector 
space of matrices of the form Xq + q, where q varies over pure quaternions. Then, the centralizer 
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of A'' in so(n,R) is the space of block diagonal matrices Bd(A5 + q, A), where q varies over pure 
quaternions and A varies over 2x2 skcw-symnictric block diagonal matrices. 
Hence, we can find a smooth g{p) e 50(n,R), and a smooth such that: 

N{p) ^ g{p)-'T{^iip))gip), with (8) 

T{^l)^BdiXi^i)q + qi^i),Ai^l)). 

The versal deformation T(fj,) is not universal (which means that ^{p) is not uniquely determined 
by N{p)), however, the nondiagonal eigenvalues of T(ju) are given by the Formula ([7]). It follows 
that q is determined modulo conjugation by a unit quaternion. On the other hand, the functions 
A(/i), A(/i) are smooth and A(/z) is nonzero. 

A. 3 Codimension of double and triple eigenvalues 

Lemma 2 We have the following: 

(i) the set oj skew symmetric matrices with a double eigenvalue is an algebraic subset of codimen- 
sion 3 in skew symmetric matrices; 

(a) the set of skew symmetric matrices with a triple eigenvalue is an algebraic subset of codimen- 
sion 8 in skew symmetric matrices. 

The proof of (i) is given in the appendix of |10] . The proof of (ii) given hereafter is a general- 
ization. We restrict ourself to the even dimensional case so(2n), the odd dimensional case being 
similar. 

We consider the set V of block-diagonal matrices D of the form 

D = Bd(Q; J, aJ, aJ, a^J, . . . , UnJ), 

of dimension N = 2n, with J ~ 0^ ) ^^'^ show only that the union of the orbits under 

orthogonal conjugation of the elements of V has codimension 8 at least. To do this we consider 
generic elements of V only: for non-generic elements the dimension of the orbit is smaller. 

To compute the dimension of the orbit Od of D it is enough to compute the dimension of the 
stabilizer G of D, and then to compute the dimension of the Lie algebra C =:Lie(G'), which is just 
the centralizer C oi D. 

By a direct computation one gets (if _D is a generic element) that elements C of C are of the 
form 

C = Bd(Ai,A), 

where is 6 x 6 and A is block diagonal with 2x2 blocks. Both Ai and A are skew-symmetric 
and Ai is of the form 

a2 J 





and 



Si J = 



-1^.3 



Then, dim(C) =dim(G) = 77,-3 + 9 = n + 6. Therefore dim(C'£)) = m — n — 6, where m = n(2n — 1) 
is the dimension of so(2n). The dimension of 2? is n — 2. Hence the dimension of the union of the 
orbits through points ofPism — n — 6-f-?i — 2 = ?7i — 8. 
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A. 4 Proof of Lemma [Tt Genericity of (R) 

Wc consider the set S of corank-2 sub-Riemannian metrics on a fixed manifold M, equipped with 
the Whitney topology. The result being essentially local, we may assume that M is an open set 
of R", with global coordinates ^, and that our sub-Riemannian metrics are globally specified by 
an orthonormal frame, i.e. s = (Fi, ...,Fp). 

For the moment, we fix s € 5. We consider two independant one forms wi,W2 on Af, that 
vanish on A, and wc set Li = dco^^, and Li is the skew-symmetric matrix defined by Li via the 
metric, and moreover we impose (as in 12. 2p that Li(^),i2(0 Hilbert-Schmidt-orthonormal. 
The matrices £1,^2 arc defined uniquely modulo a rotation Li = cos(q!(^))Li + sin(Q:(^))L2, 
L2 = — sin(a(^))Li + cos(a(^))i2- They are the same as the matrices Li in Section \^A\ and they 
meet: {L,)kj = uj,i[Fk,Fi]) ^ duj,{Fk,Fi). 

In coordinates, we set z = {^,0), and A{z) = cos(0)ii(^) -I- sin(0)i2(O- 

We fix a point a point Zq = {9o,£,q) G tig C x M, and we work in a neighborhood of zq. By 
what has just beeen said, we can perform a constant rotation to have 0q = 0. 

Local coordinates ^ ~ (x,y) in M around ^0 can be found, with xq = 0, yo = 0, such that: 
l-^'(^o) = ^,* = l,...,P 

2. WjX^o) = duj - x'Lj{£^o)dx, j = 1,2 

3. A{zo) — ii(^o) is 2 X 2 block diagonal with decrcasingly ordered (moduli of) eigenvalues. (For 
this last point, we use an (irrelevant) rotation in the distribution A^j, , i.e. a constant rotation of 
the orthonormal frame) 

Remark 4 Note that at a point zq — (^Oi^o) G Us, the two (moduli of) highest eigenvalues of 
A{zo) are equal. However, the whole construction here holds at each point of x M. 

In these coordinates, we can write (locally) s in the following form: Fi{^) = (e^ -|- B^^)-^ + 
x' Li{^o)ei-^ + x' L2{^o)ef^ + 0^{^), where = (0, .., 1, .., 0) is the i*'' coordinate vector in Rp, 
where is a px n matrix, and O^(^) is a term of order 2 in ^, i.e. O^(^) is in X^, where T is the 
ideal of smooth germs at in M", generated by the components £,i. 

This choice of notations for the vector fields Fi is adapted to the transversality arguments we 
want to apply later. Note that Li(^o), B^ are, in coordinates, components of the one-jet j^s(^o) 
of s at ^0 • 

Define the pxp matrix J7'' by U'[j = Bf ^, and by Jli , f22 the skew-symmetric matrices associated 
with the 2-forms dwi|A(Co)j c^'^2|a(Co) in the chosen coordinates. It is not hard to compute the 
tangent mappings TLi(^o) a-nd TL2{S,o) (we temporarily write TL{^q) and f2 for convenience): 

TL{^o)(er) ^W'n^n'U'', r = l,...,n. (9) 

For this, one just uses dod — 0, and TLk.i{£,o){er) = dLd{TFk{er), Fi) + dLd{Fk,TFi{er)), where 
dio stands for dtoi or du!2- 

On the other hand, we have, using the versality theorem in a neighborhood of Zq: 

A{z) = cos(6l)ii(e) + sin(6l)i2(e) = H{z)Bd{X{z)q + q{z), A{z))H'{z), (10) 

in which we already assumed that = 0, and the coordinates ^ = (x, y) were already chosen 
for A{zq) to be diagonal. Also, H{zq) = Id. 
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Remark 5 1. The decomposition iJOj) is not unique: the quaternion q{z) is defined modulo 
conjugation by a unit quaternion, q{z) = qi{z)q{z)qi{z)^^ . However, the tangent mapping Tq{z) 
is changed for Tq{z) ~ [T qi{z) , q{z)] + qi{z)Tq{z)qi{z)^^ . But on Us, q{0) ~ 0, hence the rank of 
Tq(zo) remains unchanged. 

2. The decomposition can easily be made unique, by making (following Arnold f^) some par- 
ticular choice of a (mini) transversal to the centraliser of A{zq). For instance, one could chose the 
(Hilbert- Schmidt) orthogonal supplement to the centralizer of A{zq) through A{zq). 

Let Hq : so{n) — > Q ~ M'', be the projection associating to the matrices, the quaternionic 
components of the first 4x4 block on the diagonal. 
By (dni), we have: 

Tqizo){V) = IIqTA{zo){V) + IlQ[TH'izo){V), Bd(zo)]. (11) 

We can consider the fiber mapping ttq : J^SxS^ Q x M{3,n + 1), ttq^zq ■ (j^s(Co), ^^o) ^ 
{q{zo),Tq{zo)) {Ai{3,n+ 1) being the set of 3 x [n + 1) real matrices), 

^Q,,„ (ii, L2, B\ i = 1, r) = {nQ(A(zo)), Hq o TA{zo)}. 

The following lemma is an easy consequence of (|9]), (even easier to prove if one considers thal|f| 
ri = rii = Li, is 2 X 2 block diagonal, the 2 first blocks being both nonzero): 

Lemma 3 The mapping t^q.zq is a linear submersion. 

It follows from Lemma[3]that the mapping p : J^SxS^ — > R'^ x A^(3, n + 1), (zo, Li, L2,B^,i = 
1, ...,p) — > Tq{zo) is a submersion. 

The codimension dp of the algebraic set of 3 x {n + 1) matrices that have corank 1 at least 
is do = {n — 1) [product of coranks in the 3 x (n + 1) matrices]. By Lemma ^ the set of 
skew-symmetric matrices that have double maximum eigenvalue is di = 3. Therefore, by the 
transversality theorems [5], there is a residual subset of the set of sub-Riemannian metrics, for 
which the codimension of the set of z = {9,0 in 5*^ x M where A{z) corresponds to a double 
eigenvalue, an property (R) holds at {0,£,), is a stratified set of codimension do + di = n + 2. 

A. 5 Volume of the unit ball 

We keep the notations of Section [3] 

£t{p'o.Pl)=St{p'o:r,e)^£i{tpltr,e), 
ant tcut = A^{9)/r. The domain of £t for the unit ball is 

Ds, - {{p^o^r,0,t) I e£ [0,2tt], e SP-\t e mm{l,t,ut), re [0,+^[}, 

where Sp~^ denotes the unit Euclidean sphere in M.P and later B denotes the unit Euclidean ball 
W. In this formula, the boundary of this set in the variables r, t is parametrized by r. Equivalcntly 
if we parametrize this boundary by t we get, 

Ds^={{p-o,r,0,t) I 0e [0,27r], G t e [0, 1], r^A^{d)/t}. 

Now set p = tpo, f — tr. This implies for the domain Dg-^ oi £i{p,r,9), 

*Note that in fact, in the chosen coordinates, Qi = Li since Fi{za) = on Azq 



11 



D£,={{p,p,e) I 0e [0,2^], pes, f e[o,A^{e)]} . 

The volume of the unit ball of the nilpotcnt approximation at ^ is 



V^= I Popp = / A dy, 



that is 

= / / / Jsiip.f. 
Jo Jo Jb 

where J^^ is the Jacobian determinant of £i{p,r,0). 



e)dpdf de, 
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